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Abstract 

The nonperturbative contribution, suppressed by powers of the charm quark mass 
rric, in the inclusive decay B Xs'^ is analyzed in the context of heavy quark 
effective field theory (HQEFT). According to previous analyses, the leading effects 
of 0(l/m^) arise from an external gluon attached to a charm quark loop and can 
be expressed in terms of the chromomagnetic interaction of the h quark inside the 
B meson. This is also true at leading order in the HQEFT approach. However the 
structure of higher-dimensional operators is different because the effects of external 
gluons alone do not give a complete set of operators. A systematic method to derive 
all the high-dimensional operators can be obtained in the HQEFT scheme using the 
operator product expansion. 

Key words: nonperturbative effects, heavy quark effective theory, operator 
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The inclusive radiative decay B Xs7 has been considerably investigated 
experimentally and theoretically. The CLEO collaboration [1] reported the 
branching ratio for this decay to be B{B X,7) = (2.32±0.57±0.35) x 10"^. 
Theoretically this type of rare decay modes can be a sensitive probe for new 
physics beyond the standard model. And the photon energy spectrum is also 
interesting in understanding CP violation and new physics in B decays [2]. In 
order to probe new physics, there should be a precise theoretical calculation 
of the branching ratio in the standard model. For the decay B — >■ Xs7, the full 
next-to-leading order calculation in the standard model is completed combin- 
ing matching conditions [3], matrix elements [4] and anomalous dimensions [5] 
in the effective Hamiltonian approach. 
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Since the advent of the heavy quark effective field theory (HQEFT), the un- 
derstanding of the decay of B mesons has come to a more advanced level. Ap- 
plying the HQEFT to B mesons, the inclusive decay rates such as B ^ Xcii' 
and B —>■ Xgj can be obtained by simply calculating the corresponding quark 
decay rates at leading order [6]. The corrections to the leading result can be 
systematically incorporated in the scheme of the HQEFT. These corrections 
consist of a double series expansion both in the strong coupling constant ag 
and in the inverse power of the b quark mass, mf,. Systematic analyses of 
the nonperturbative effects suppressed by powers of rrif, have been considered 
extensively in various inclusive B decays [7]. With all these ingredients the 
theoretical uncertainty is reduced to a 10% level. The object of this letter is 
to elaborate on the contribution of high- dimensional operators, especially of 
order 0(l/m^) in the context of the HQEFT. 

The theoretical framework to analyze the inclusive decay B Xg^ is to use 
the effective weak Hamiltonian [8] which is given by 

AG ^ 



where Oj are various operators. The operator O2 is the four-quark opera- 
tor given by O2 — {sLaliJ>Lp)icLpl^CLa) in which a, (3 denote color indices. 

The operator Oi differs from O2 only in the way the color indices are con- 
tracted. Operators O3 to Og are four-quark operators including all flavors 
below the renormalization scale ^ of order tth,. The operator O7 is given as 
O7 = {e/167i'^)mbSLcr^'^F^i,bR, and Og is obtained from O7 by replacing eF^j, 
by gG^^,. 

The contribution of these operators to the inclusive B Xgj decay rate is sys- 
tematically calculable. In the leading log approximation, the matrix element of 
C'j{ijl)0'j{ii) dominates in evaluating the decay rate B — > Xg^ for large enough 
photon energies. At next-to-leading order, other operators can contribute. In 
analyzing the contribution of other operators, when it is possible to use the 
operator product expansion (OPE), higher-order processes can be described 
by products of some coefficients which are calculable in perturbation theory, 
and local operators whose matrix elements represent long-distance effects. 

When operators in iJcfr other than O7 are included, if we try to calculate 
perturbatively the decay rate for B — > Xs7, many problems arise because the 
decay rate receives contribution from those in which the photon couples to 
light quarks. Actually the amplitude for h sgj with a light quark loop is 
nonlocal. Therefore for light quarks coupled to a photon, the idea of the OPE 
is not applicable. That is, there is no OPE to parameterize nonperturbative 
effects from the photon coupling to light quarks in terms of B meson matrix 
elements of local operators. 
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However there are nonperturbative effects arising from the photon couphng 
to the charm quarlc, and it is possible to have B meson matrix elements of 
local operators which are suppressed by (Aqcd/'^T'c)^ rather than (Aqcd/''716)^. 
We can calculate the nonperturbative effects suppressed by I /ml because 
the charm quark is heavy {rric ^ Aqcd)- It is exactly the reason why we 
can employ the HQEFT idea in order to obtain high-dimensional operators 
systematically. 

Naively we expect that the high-dimensional operators, which result in non- 
perturbative effects, are suppressed by powers of rrib since the only physically 
relevant scale in the inclusive decay B — > Xs'~f is the mass of the b quark. 
However Voloshin [9] has observed that there arc also nonperturbative effects 
which are suppressed by powers of rric, the charm quark mass. This nonpertur- 
bative effect arises from the operator O2 as shown in Fig. 1, which is produced 
by the quantum fluctuation of the cc pair emitting a hard photon, and an ex- 
ternal soft gluon is attached to the charm quark loop. Here the blob represents 
the operator O2. 




Fig. 1. Feynman diagram for b sg'j from O2 in which a photon and a gluon are 
attached to the c quark loop. The diagram with a photon and a gluon switched is 
omitted. 

If it is possible to use the OPE in obtaining high-dimensional operators from 
the effective Hamiltonian in Eq. (1), the contribution from the operator O2 is 
the most important one because the coefficient C2 is large compared to C-^-Cq. 
The high-dimensional operators arising from the operator Oi starts from the 
operator with more than one gluon due to the color structure, hence more 
suppressed compared to the operators obtained from O2. 

To first order in the gluon momentum, the new effective Hamiltonian for 
b — >• sjg is written as 

Here 0^/3 = C^^Ta is the gluon field strength tensor, F^x is the electromagnetic 
field strength tensor, and = d"' + ieQcA"' + igG°' is the covariant derivative. 
Qc = 2/3 is the electric charge of the charm quark. And the sign convention is 
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g0i23 _ _|_]^ There has been much effort to estimate the size of this new type of 
nonperturbative effects in various B decays such as B ^ ^sl, B Xgi'^i' 
[10-13]. 

If we keep more powers of external gluon momentum, or if we attach more 
giuons to the charm quark loop, we can obtain higher-dimensional operators. 
These higher-dimensionaf operators are suppressed by more powers of rric com- 
pared to the leading operator given in Eq. (2). These higher-order effects are 
quite small compared to the contribution of the leading operator, not because 
the expansion parameter mi,h.qcY>/m1 is small, but because the numerical co- 
efficients of these operators are small [10,11]. 

Now we consider how we can implement the idea of the HQEFT and the 
OPE to derive high-dimensional operators suppressed by powers of rric in a 
systematic way. In the heavy quark limit m;, — > oo, the h quark interacts with 
the light degrees of freedom, which is characterized by the residual momentum 
k of the h quark (p^^ = m^v^ + k^). As the operator O2 acts, the b quark turns 
into an s quark and a cc pair. When the virtual cc pair emits a hard photon 
and soft giuons, this process induces high-dimensional operators. Now let us 
also take the heavy quark limit mc —> 00 while rric/ mi, held fixed. Then when 
the b quark turns into a cc pair, the cc pair is still immersed in the same light 
degrees of freedom from the b quark. Therefore the c quark also carries the 
residual momentum k of the original b quark. The fact that the charm quark 
should have a residual momentum is manifest if we consider the corresponding 
Feynman diagram in the full theory with W exchange. 

From this process mentioned above, we can systematically obtain all the high- 
dimensional operators in powers of l/rric by expanding the residual momentum 
in the c quark and using the OPE. There are also other contributions such as 
the effect of the residual momentum in the s quark, but the operators obtained 
from the residual momentum of the s quark are suppressed by powers of m^. 
Since we are interested in the operators suppressed by powers of nic here, we 
will neglect all the operators suppressed by 171^. And at this leading order, the 
b field can be regarded as a heavy field by with momentum = rrihV^ + k^ in 
the HQEFT. 

These high-dimensional operators include the operators obtained in Ref. [9- 
13], which correspond to the so-called "one-gluon matrix element" [7]. The 
operators from one-gluon matrix element are those in which an external gluon 
is attached to the charm quark loop. In order to make the resultant operator 
gauge invariant, we combine one gluon momentum operator and the exter- 
nal gluon field to form the gluon field strength tensor Gajs- Therefore these 
operators are proportional to Gap, and its derivatives if there are additional 
gluon momenta at higher orders. If we attach more external giuons, we obtain 
high-dimensional operators in powers of Gap. However, in the HQEFT, there 
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are also other operators which arc not proportional to Gap and its deriva- 
tives, and these operators can be obtained systematically using the HQEFT 
and the OPE. Interestingly enough, the leading contribution in the HQEFT 
is the same as the operator from the one-gluon matrix element. Therefore the 
numerical analysis for the contribution of the leading operator at 0(l/m^) 
in Ref. [9-11] does not change. The contribution of the leading operator in 
Eq. (2) compared to the contribution of 0^ alone is about 5V /V w 3% [13]. 

In order to systematically calculate high-dimensional operators suppressed by 
powers of iric in the HQEFT, consider the Feynman diagram shown in Fig. 2. 
Though the operator O2 is shown in Fig. 2, note that the operator Oi should 
be included also when there are two or more soft gluons are emitted. However 
the contribution of Oi appears at next-leading order compared to that from 
02- The integral in evaluating the Feynman diagram in Fig. 2 can be written 
as 

J (27r)^ '-f — § + f — rric f + f — rric 



where q denotes the photon momentum. 




Fig. 2. Feynman diagram for 6 ^ 57 to obtain high-dimensional operators from the 
expansion of the residual momentum in the charm quark loop using the OPE. 

Now we expand /'^^ in powers of the residual momentum k and collect the 
terms which are suppressed by powers of rric- We use the operator product 
expansion in order to obtain the corresponding operators. That is, we replace 

k'^ by H^ = iD'^, where D^^ is the covariant derivative acting on the b field. 
The operator H^ satisfies the commutation relation [n^,n'^] = —igG^^. The 
amplitude for 6 — > 57 -|- many gluons from Fig. 2 can be written in terms 
of I^" as 

iM' = ^V:,V^C2{sl'^^,{l - 75)&). (4) 



At a first glance, the integral in Eq. (3) seems to be made independent of the 
residual momentum when we shift the loop momentum variable Ihy l+k, hence 
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appear no high-dimensional operators. But this shift of the integration variable 
is not allowed in the integral because the integral as such is quadratically 
divergent. However, if we employ the dimensional regularization to regulate 
the integral, only the 75-dependent part of the integral should be treated 
carefully. 

The 75-independent integral can be shifted to be independent of the residual 
momentum, hence producing no high- dimensional operators. Furthermore af- 
ter the shift of the integration variable, the 75-independent integral vanishes 
for the on-shell photon. Therefore to all orders there are no high-dimensional 
operators from the 75-independent part. It means that high-dimensional op- 
erators suppressed by powers of rric in B ^ Xs'~f should always include the 
Levi-Civita e tensor which appears by taking the trace of Dirac matrices in- 
volving 75 in Eq. (3). 

In order to calculate the remaining 75-dependent part which is suppressed by 
powers of rric, we first expand the integrand in powers of k, and we can shift 
the integration variable only when the integral is guaranteed to be convergent. 
In the resultant expression, we get high-dimensional operators by replacing k 
with the covariant derivative iD. 

With this in mind, let us consider 7^'^ in Eq. (3). When we first replace k by 
the covariant derivative 11 = iD, and expand 7^^ in powers of 11, we get 



K"^^) m,n=0 



, n+l 



X tr 



g-^+yL + m,) (2n . (/ - g) + - ^G^^a'^'^y 
X 7^(/ + + me) (2n • Z + - ^Gpra'^fril - Ts) 



(5) 



where o""^ = i['~f°', 7^]/2. Here the appearance of GajS comes from thepl/t term 
when we expand the denominator in Eq. (3) with the identity 



(6) 



From Eq. (5) we can systematically obtain high-dimensional operators by a 
series expansion with respect to m and n. Here I replace the residual mo- 
mentum by the covariant derivative first and treat the order of the products 
carefully as in Eq. (6). Because of this procedure, there appears Gap in Eq. (5). 
In the literature [7], there is another way of calculating the integral. First we 
expand the diagram in powers of the residual momentum k, regarding it as a 
number, i.e., 1^ = H^. Then there appears no G^is in the expansion. In order 
to include the effect from the gluon field strength tensor Gap, we attach an 
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external gluon and extract the term proportional to Gajj- It is called a "one- 
gluon matrix element". However it turns out that these two approaches are 
equivalent. 

It is complicated to obtain a simple and general form for the series expansion 
in Eq. (5). However we can calculate the first few high-dimensional opera- 
tors explicitly in order to see the structure of them. The zcroth-ordcr term 
(m = n = 0) does not contribute to the high-dimensional operators suppressed 
by rric- Also there is no term linear in H. At second order in H, there are no 
symmetric terms quadratic in H such as H"H^-|-H^n" or H^. The only surviv- 
ing terms up to quadratic terms in H are those with Gas- These terms come 
from the series expansion in Eq. (5) with m -\- n — 1 and 2 respectively. The 
result is written as 



egQc 
247r2m? 



,>^-^<^q^qPG^p. (7) 



Note that all the quadratic contributions which do not involve Gaj} explicitly 
in Eq. (5) cancel out, leaving only the one-gluon matrix element. The reason 
why the term with Gap alone survives is clear when we look at the diagram in 
Fig. 2. If a single gluon is attached to the loop, the loop has the structure of 
a vector-vector-axial triangle and only the term with the Levi-Civita e tensor 
and survives. The effective operator from Eq. (7) can be derived in a 
straightforward way, and is given in Eq. (2). Therefore the leading term of 
order coincides with the result by Voloshin [9], which is obtained by 

attaching a soft gluon to the c quark loop. 

Now let us consider higher-order corrections to the leading term given in 
Eq. (7). The operators of 0(l/m^) are obtained from the terms with m-\-n — 2 
and 3 respectively in Eq. (5). The result is written as 

egQc -Aa..^^^/3[^.n,G^^]. (8) 



ISOTT^m^ 



Recall that the covariant derivative operator H is applied to the h field. Then 
the commutator in Eq. (8) can be written as 

[q-Yl,G^p]h^{q-YlG^p)h, (9) 



where the covariant derivative H in the right-hand side only applies to G^p- 
Of course, it should be understood that this holds only for the derivative in 
the covariant derivative H. Using the identity Eq. (9), Eq. (8) can be written 
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as 



egQc 



(10) 



ISOTT^m^ 



where 11 applies only to the gluon field strength tensor. This is where the 
HQEFT result is different from the result by attaching soft gluons to the charm 
quark loop. To be more quantitative, the HQEFT result is half the higher- 
dimensional contribution derived from the one-gluon matrix element and its 
derivatives, which makes the contribution of higher-dimensional operators of 
0{l/mf) smaller than estimated before [10,11]. 

There are also other operators which cannot be obtained by attaching ex- 
ternal gluons. These operators can be obtained in the HQEFT scheme. For 
example, there are higher-dimensional operators of order H^/m^, though their 
contributions are suppressed by A/mj compared to the operators of order 
iP/m^. These terms arise from m + n = 1, 2, 3 respectively. After a long but 
straightforward calculation, the result is given by 



In evaluating the contribution of higher- dimensional operators of order O(H^), 
Eqs. (8) and (11), to the decay rate, there is one more caveat that we have 
to include the contribution of both the operators Oi and O2. For the leading 
term, given in Eq. (7), with one gluon field strength tensor, only the operator 
O2 contributes since it has the correct color structure. However, as mentioned 
before, when we include the emission of two or more gluons, the contribution 
from the operator Oi should also be considered. We have to separate these 
higher- dimensional operators into the color-octet part for O2 and the color- 
singlet part for Oi. But since these contributions are numerically insignificant, 
it will not be pursued here any more. 

Implementing the HQEFT idea to the charm quark which is immersed in the 
same light degrees of freedom of the decaying b quark, we can systematically 
obtain a complete set of high-dimensional operators suppressed by powers of 
rric. The effect of attaching external gluons gives the correct leading term, 
but it does not include all the higher-dimensional operators. In other words, 
attaching external gluons corresponds to the expansion in powers of the gluon 
field strength tensor Ga/s and its derivatives. The resultant operators for n 
attached gluons are symbolically of the form with some derivatives. In 




(11) 
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order to obtain a complete set of high-dimensional operators, the HQEFT idea 
can be applied to the decay B — > Xgj. 

The numerical estimate of the higher-dimensional operators are negligible 
compared to the leading operator in Eq. (2), which already contributes about 
3% to the decay rate. As pointed out in Refs. [10,11,13], it is not because 
the expansion parameter nibAQCB / 'itT'I is small, but because the coefficients 
are small. Even though higher-dimensional operators are numerically insignif- 
icant, the idea of the HQEFT offers a systematic way of obtaining a complete 
set of high- dimensional operators. 

Conceptually there is an advantage in considering the nonperturbative con- 
tributions using the HQEFT idea. The effect of the operators suppressed by 
powers of rric is nonperturbative. The reason why we could systematically clas- 
sify high-dimensional operators is because the charm quark is heavy compared 
to the QCD scale Aqcd- If the charm quark is not heavy, it is impossible to 
analyze the l/rric effects using the OPE, and this effect belongs to a purely 
nonperturbative regime. And it is not surprising to have nonperturbative con- 
tributions governed by the scale rric as well as rrih in B decays since these 
effects arise simply because the c quark is also heavy. In the same line of rea- 
soning, it is impossible to analyze the contribution of the up quark loop from 
a similar process in inverse powers of the up quark mass using the OPE since 
this contribution is genuinely nonperturbative. 
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